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ABSTRACT 

In this paper we analyse and compare the clustering of young stars in Chamaeleon I 
and Taurus. We compute the mean surface-density of companion stars ./V as a function 
of angular displacement 9 from each star. We then fit N(9) with two simultaneous 
power laws, i.e. N(9) ~ K bin 6>~ /3bin + A" clu 6> _/3clu . For Chamaeleon I, we obtain (3 bin — 
1.97 ± 0.07 and /? clu = 0.28 ± 0.06, with the elbow at elb ~ 0.011° ± 0.004°. For 
Taurus, we obtain bin = 2.02 ± 0.04 and f3 clu = 0.87 ± 0.01, with the elbow at elh ~ 
0.013° ± 0.003°. For both star clusters the observational data make quite large (~ 5a) 
systematic excursions from the best fitting curve in the binary regime (9 < 9 elb ). 
These excursions are visible also in the data used by Larson and Simon, and may be 
attributable to evolutionary effects of the types discussed recently by Nakajima et al. 
and Bate et al. In the clustering regime (9 > 9 elb ) the data conform to the best fitting 
curve very well, but the /3 clu values we obtain differ significantly from those obtained 
by other workers. These differences are due partly to the use of different samples, and 
partly to different methods of analysis. 

We also calculate the box-dimensions for the two star clusters: for Chamaeleon I 
we obtain D box ~ 1.51 ±0.12, and for Taurus D box ~ 1.39 ±0.01. However, the limited 
dynamic range makes these estimates simply descriptors of the large-scale clustering, 
and not admissible evidence for fractality. 

We propose two algorithms for objectively generating maps of constant stellar 
surface-density in young star clusters. Such maps are useful for comparison with 
molecular-line and dust-continuum maps of star-forming clouds, and with the results 
of numerical simulations of star formation. They are also useful because they retain 
information which is suppressed in the evaluation of N(9). Algorithm I (SCATTER) 
uses a universal smoothing length, and therefore has a restricted dynamic range, but 
it is implicitly normalized. Algorithm II (GATHER) uses a local smoothing length, 
which gives it much greater dynamic range, but it has to be normalized explicitly. 
Both algorithms appear to capture well the features which the human eye sees. We 
are exploring ways of analyzing such maps to discriminate between fractal structure 
and single-level clustering, and to determine the degree of central condensation in 
small- N clusters. 

Key words: binaries: general - stars: formation - ISM: clouds 



1 INTRODUCTION 
1.1 Motivation 

As the number of young stars detected in nearby star- 
formation regions (SFRs) has increased over the last decade, 
it has become possible to study the clustering properties in 
these regions flGomez et al. 1993) ; [Larson 1995] | imon 19971 



Nakajima et al. 1998; Bate, Clarke & McCaughrean 199£ 



* e-mail: philip.gladwin@astro.cf.ac.uk 

t e-mail: spyridon.kitsionas@astro.cf.ac.uk 

t e-mail: henri.boffin@astro.cf.ac.uk 

§ e-mail: ant.whitworth@astro.cf.ac.uk 



Brandner Kohler 1998| ). The manner in which young stars 
cluster places important constraints on theories of star for- 
mation. Numerical simulations of star formation are now 
approaching a level of sophistication which enables them to 
make detailed predictions about the density- and velocity- 
fields in star-forming gas and the distributions of the result- 
ing stars (e.g. Turner et al. (1995), Whitworth et al. (1995), 
Bonnell et al. (1997), Klessen (1997), Bonnell, Bate & Zin- 
necker (1998), Boffin et al. (1998)). Since there is a great 
diversity of star-forming molecular clouds and young star 
clusters, useful comparisons between observation and sim- 
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ulation can probably only be made using global statistical 
descriptors, and this paper seeks to develop such descriptors. 

1.2 Background 

Speckle interferometry, lunar occulation and spectroscopic 
surveys have enabled astronomers to estimate the binary 
fraction in young stellar populations. In general, the pro- 
portion of binaries is at least as high as for mature stars in 
the field, and in Taurus it appears to be even higher. At the 
same time, near infrared imaging and X-ray surveys of SFRs 
have revealed many additional young stars which were not 
detected by the earlier emission-line searches. 

Extending earlier work by Gomez et al. (1993), Larson 
(1995) has computed the mean surface-density of compan- 
ions per star, N(8), as a function of angular separation, 
8, in Taurus. Since two distinct power laws of the form 
N(8) ~ 6~P are necessary to fit his data, Larson concludes 
that there are two different clustering regimes. The smaller 
separations define the regime in which binary and higher- 
order multiple systems are found; in this regime he finds 
/3 b in ~ 2.15. The larger separations define the regime in 
which stars appear to be hierarchically clustered, and here 
he finds /? olu ~ 0.62. Larson notes that the break between 
the two regimes (hereafter referred to as the elbow) occurs 
at s c]b ~ 0.04 pc, and suggests that this is the Jeans length 
in the cloud cores from which the stars have formed. 

Simon (1997) has repeated this analysis for Taurus, and 
extended it to the Orion- Trapezium and p Ophiuchus SFRs. 
He finds values of /3 b i„ and /3 c i u similar to those found by Lar- 
son (1995). However, the elbow shifts to smaller separations 
going from Taurus (s o]b — 0.06 pc), through p Ophiuchus 
(0.024 pc) to Orion- Trapezium (0.002 pc). Simon argues 
that it is difficult to explain this shift simply on the basis of 
a 30-fold variation in the Jeans length. He points out that 
in most cases the elbow is determined largely by the back- 
ground surface-density of the overall cluster; this is reflected 
in the fact that s c i b correlates with the mean separation be- 
tween stars. Essentially the same conclusion was reached by 
Kitsionas, Gladwin & Whit worth (1998). 

Nakajima et al. (1998) have calculated N(9) for the 
Orion, Ophiuchus, Chamaeleon, Vela and Lupus SFRs. For 
SFRs where two power laws are required to fit the data, they 
obtain f3 hin ~ 2.0 and 0.1 <, P a i u <, 0.8. They show that the 
elbow occurs at a separation roughly equal to one tenth of 
the mean nearest-neighbour distance. They also point out 
that structure which is present when the stars first form 
will subsequently be erased by random migration due to the 
velocity dispersion, which is typically ~ 1 km s _1 . Thus, 
as unbound clusters disperse, the elbow should shift to large 
separations. In addition, some SFRs show evidence that star 
formation has gone on for several Myrs; and so there may 
be a young, highly clustered population mixed in with an 
older, less highly clustered population. 



Bate, Clarke & McCaughrean (1998) have recently un 



dertaken a comprehensive evaluation of the factors influ- 
encing N(8) and its interpretation. They show that N(8) 
can be strongly affected by boundary effects. They point 
out that even if N(8) can be fit with a non-integer power- 
law, this does not necessarily imply fractality. They illus- 
trate this by considering non-fractal distributions for which 
N(9) can nonetheless be fit with a non-integer power-law. In 



particular, a simple distribution of single-level clusters in a 
lower-density background can produce an N(8) which is well 
fit by a non-integer power-law. They also demonstrate that 
evolutionary effects - dispersal of unbound clusters due to 
their intrinsic velocity dispersion, evaporation of stars from 
small- N bound clusters, tidal dissolution of wide binaries - 
should all act to suppress N(8) in the vicinity of the elbow. 

Whitworth, Boffin & Francis (1998) have argued that 
binary formation via dynamical fragmentation only occurs 
when the gas becomes thermally coupled to the dust. This 
leads to a critical density below which binaries are unlikely 
to form via dynamical fragmentation, and hence to a maxi- 
mum initial binary separation ~ 0.05 pc. They suggest that 
this may contribute to the break in N(8) at around this sep- 
aration in sparse star clusters like Taurus, where the overall 
density is too low for chance projections to be important. 



1.3 Plan 

In this paper we concentrate on the Chamaeleon I SFR, a 
regi on of low mass star formation at a distance of ~ 140 
pc (Schwartz 1992). In Section we evaluate N(6) and the 



box-dimension of the star field, and in Section 3 we review 
briefly the procedure used by Gomez et al. (1993) to generate 
surface-density maps of young star clusters. In Sections 4 
and 5 we describe and evaluate two new algorithms that 
create contour maps of stellar surface-density. In Section 6 
we compare our results for Chamaeleon I with Taurus (also 
at a distance of ~ 140 pc (Elias 1978)), and in Section 7 we 
summarize our conclusions. 

For Chamaeleon I we use the positions of 137 optically 
visible T-Tauri stars compiled by Lawson, Feigelson and 
Huenemoerder (1996) and K-Band observations of binary 
systems by Ghez et al. (1997). For Taurus, positions of 216 
optically visible pre-main-sequence stars have been compiled 
from the sources listed in Table ^ these are essentially the 
same data used by Gomez et al. (1993) but supplemented 
by the binary searches of Ghez et al. (1993), Leinert et al. 
(1993) and Simon et al. (1995). 



2 SURFACE-DENSITY OF COMPANIONS 

In order to study the distribution of young stars in Taurus, 
Gomez et al. (1993) evaluated the two point angular corre- 
lation function. Larson (1995) instead evaluated the mean 
surface-density of companions per star, N(6), as a function 
of angular separation, 8; this measure avoids the need to de- 
fine a finite background surface-density. Simon (1997) used 
both methods to study the distribution of stars in three 
different SFRs, and in each case calculated N(8) out to sep- 
arations ~ jq of the overall angular extent of the region, 
where the two methods are in good agreement. Nakajima et 
al. (1998) also limited their fits to 8 < 6> max with max less 
than the overall size of the SFR being analyzed; but they 
do not divulge how (9 max was chosen, and our numerical ex- 
periments suggest that the fit can be quite sensitive to the 
choice of 9 m „. 
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Table 1. Summary of references used to compile the Taurus catalogue. 



Reference Number of stars Survey Type Notes 



Herbig and Bell (1988) 


118 


Slit spectroscophy 


(1) 


Walter et al. (1988) 


14 


X-ray survey 


(1) 


Hartmann et al. (1991) 


10 


Proper motion survey 


(1),(2) 


Gomez et al. (1992) 


2 


Proper motion survey 


(1) 


Briceno et al. (1993) 


12 


Ha prism survey 


(1) 


Leincrt et al. (1993) 


167 


nIR speckle interferomctry 




Ghez et al. (1993) 


67 


Speckle imaging survey 




Simon et al. (1995) 


77 


IR, lunar occultation survey 





Notes: (1) as described in §2 of Gomez et al. (1993), (2) includes all the stars of 
their Table 1. 



2.1 Data reduction 

To compute N(9), we take each star n in turn, and divide 
the surrounding area of sky into a set of annuli, by drawing 
circles of radius 9i centered on star n, with 9i = 2#i_i and 
i > 1. 9o is chosen to be well below the smallest separation 
in the sample. Next we count the number of companion stars 
Af n ,i in each annulus. Then, 



N($4 



-i) 



(1) 



where 9i = (9i + /2 and 7V* t ot is the total number of 
stars in the cluster. 

Because of discreteness and small-number statistics, the 
results are quite sensitive to the choice of 6q. Therefore we 
repeat this procedure twice, each time increasing 9o by 2 1 / 3 . 
This increases the number of plotted points threefold. It 
appears to resolve all the significant divergences between 
the plotted points and the fitting curves (see, for example, 
the dip around 9 ~ 0.02° in the plotted points for Taurus 
on Figure 1), and therefore there is nothing to be gained by 
increasing the number of points any further. We believe that 
it is important to resolve these divergences, at the very least 
because they provide a graphic measure of the limitations 
of the fitting curves. It may even be the case that these 
excursions constrain the previous evolutionary history of the 
SFR (cf. Nakajima et al. 1998, Bate et al. 1998). 



2.2 Corrections for incompleteness 

Since the stars in Chamaeleon I and Taurus have not all 
been surveyed for multiplicity, and since some of those which 
have been surveyed have been surveyed to different limits 
and/or by different techniques, we have attempted to correct 
for incompleteness. To do this we assume that multiplicity 
amongst young stars is scale-invariant, in the sense that a 
faint system is as likely to be a binary as a bright one. Hence, 
in Chamaeleon I we presume that the stars not surveyed for 
companions by Ghez et al. (1997) have - on average - the 
same multiplicity as those that were surveyed. In Taurus we 
are guided by the multiplicity statistics for stars surveyed 
by Ghez et al. (1993), Leinert et al. (1993) and Simon et 
al. (1995). Systems are categorized according as they were 
surveyed in all three programmes, in two, in just one, or in 
none. They are then compensated for potential companions 
that would have been missed, on the assumption that the 
probability for companions in the ranges not surveyed is 



the same as for stars that were surveyed in these ranges. We 
emphasize that these corrections only change the coefficients 
of the fitting curves (Ky, in and K c i n in Eqn. (^)); they do not 
significantly affect the exponents (/3 b in and /3 c iu in Eqn. (^)). 

Bate et al. (1998) discuss a Monte Carlo technique to 
correct for companions beyond the boundary of the sur- 
vey region. However, this correction is only appropriate in 
cases where the cluster continues - at a comparable surface- 
density - outside the area of the survey. In the cases dis- 
cussed here this is not the case, and so no such correction is 
attempted. 



2.3 Fitting procedure 

In order to fit N(9), we follow Simon (1997) in ignoring sep- 
arations greater than one tenth of the angular extent of the 
region under consideration. We also ignore small separations 
where the counts are too low to afford reliable statistics. This 
leaves the ranges 0.000015° < 9 < 0.15° for Chamaeleon I; 
and 0.00006° < 9 < 2.2° for Taurus. The points outside 
these ranges are represented by open circles on Figure 1. 

The procedure adopted by other workers has been to fit 
the observational data piecewise with two power laws. How- 
ever, our experience has shown that the results are then 
rather sensitive to the breakpoint where the two power- 
laws meet, and there is no obvious strategy for choosing 
this breakpoint objectively. Moreover, we would - in gen- 
eral - expect there to be binaries with separations greater 
than the elbow, and chance projections at separations less 
than the elbow. Therefore our strategy has been to fit the 
observational data with two simultaneous power laws: 



N{9) = K u , 



-Phi, 



+ K cl »9 



(2) 



The elbow 9 Qih is defined as the separation at which the 
two terms on the righthand side of Eqn. (0) are equal, i.e. 



oib = exp 



In [K biu ] - In [tf clu ] 



/3ci 



(3) 



In order to identify the best fit we evaluate for each 
trial fit the probability that the observational data could be 
generated from that fit by Gaussian statistics. The best fit is 
then the one with the highest probability, and the uncertain- 
ties embrace la departures from the best fit. The data that 
we fit are the corrected data, but it should be emphasized 
that we estimate the standard deviation for the contents of 
each annulus (and hence the differential probability) using 
the raw counts, and not the corrected counts. 
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2.4 Results 

N(9) is plotted for Chamaeleon I and Taurus in Figure |l| 
The corrected points are the crosses with la error bars. For 
Chamaeleon I (lefthand panel of Figure |l]), the best fit has 
exponents f3 bill = 1.97 ±0.07 (in the binary /multiple regime 
at small separations) and /3 c i„ = 0.28 ± 0.06 (in the hier- 
archical clustering regime at large separations); the elbow 
is at 9 rAb ~ 0.011° ± 0.004°, corresponding to a projected 
separation of s c i b ~ 0.027 pc (5500 AU) at the distance of 
Chamaeleon I. For Taurus (righthand panel of Figure Q) the 
exponents are (3 bia = 2.02 ± 0.04 and /3 olu = 0.87 ± 0.01, 
with the elbow at 9 clb ~ 0.013° ± 0.003°, corresponding to 
s c ib ~ 0.032 pc (6500 AU) at the distance of Taurus. Table 
^| summarizes the values of /3 bin , /3 c i u and s i b obtained by 
Larson (1995), Simon (1997), Nakajima et al. (1998), Bate 
et al. (1998), and this work. 

We reiterate that the corrections we make for incom- 
pleteness have no significant effect on the exponents f3 bin 
and /3 c i u , or on the value of olb . They tend to increase the 
coefficients K bill and K cla by 40-50%. To demonstrate this 
we have plotted the uncorrected points as filled circles on 
Figure 1. Equally, if we count the multiple systems found 
in the surveys of Chez et al. (1993), Leinert et al. (1993) 
and Simon et al. (1995) as single stars (or "systems") for 
the purpose of computing large separations, this too has no 
significant effect on /3 clu . 

We note that in the binary regime (9 < 9 Bib ) (a) the ob- 
servational data make rather wide but apparently systematic 
excursions about the best fit line (up to 5<r); and (b) these 
same excursions are also present in the data presented by 
Larson (1995) and Simon (1997). They are more apparent 
in our plots simply because we have displayed more points. 
Such excursions might be evidence for evolutionary effects 
of the type described by Nakajima et al. (1998) and Bate 
et al. (1998), for example dispersal of unbound clusters, but 
equally they could have been present when the stars first 
formed, or they could be due to some unidentified selection 
effect. 

Duquennoy & Mayor (1991) have measured the orbital 
elements of a large field sample of mature binary systems 
having F- and G-type primaries. They show that the distri- 
bution of periods is well fitted by a Gaussian, 



Mot 

(27T) 1 /2 (J( 



exp 



-log^o [P/P a _ 



dlog 10 [P], (4) 



with P ~ 180 years and a a ~ 2.3. If we assume that all these 
systems comprise a O.8M primary and a 0.2M Q secondary 
in a circular orbit on the plane of the sky (i.e. we neglect the 
distributions of mass, eccentricity and orbital inclination), 
then the Duquennoy & Mayor distribution corresponds to 



JV DM (0) 



3/c, 



2(27r) 3 / 2 log e [10]cr o (9 2 



exp 



-91og 10 [fl/flo] 

8<7 2 



(5) 



where / comp ~ 0.9 is the mean number of companions per 
star, and 9 a ~ 0.00006°. For comparison, this curve is plot- 
ted as a thick solid line on Figures la and lb. The Gaussian 
part of Eqn. (§) has FWHM Alog lo [0] ~ 3.6, i.e. it is very 
flat. Therefore the dominant variation comes from the 9~ 2 



term, and the Gaussian part simply superimposes a subtle 
negative curvature. 

It is well known that the Duquennoy & Mayor field bi- 
nary distribution predicts fewer wide binary systems (9 <; 
0.0001°) than are observed in Taurus, and Figure lb con- 
firms this. As pointed out by Padgett, Strom & Ghez (1997), 
there are at least three evolutionary processes which might 
reduce the number of wide systems in a mature stellar pop- 
ulation, (i) Tidal disruption should preferentially reduce the 
number of wide systems, (ii) Interactions with other stars 
may also harden some systems, moving them to smaller sep- 
arations, (iii) Low-mass (substellar) companions may fade to 
invisibility. 

Additionally, it may be that different SFRs generate 
different binary distributions, and then these are all con- 
volved together to produce the Duquennoy & Mayor field 
binary distribution. There is conflicting evidence on this is- 
sue. Padgett et al. (1997) find that in the range 138 to 1050 
AU, the binary frequency in high-density SFRs in Orion is 
essentially the same as in the low-density SFRs studied by 
Reipurth & Zinnecker (1993) - and significantly higher than 
in the field. In contrast, Brandner & Kohler (1998) show that 
in Sco-Cen there are significant variations in the distribu- 
tion of binary separations in the range 15 to 450 AU, in the 
sense that sub-regions where there are more massive stars 
also tend to have closer binaries. It will be very important 
to establish how universal this trend is, for example by in- 
vestigating whether the lack of massive stars in Chamaeleon 
I and Taurus is accompanied by a relative paucity of close 
binary systems. However, the small numbers of close sys- 
tems in the data sets for Chamaeleon I and Taurus make it 
impossible to draw such a conclusion at this time. 

2.5 Projection effects 

We concur with Simon (1997) that the apparent shift m s c i b 
going from Taurus (~ 0.03 to 0.06 pc) through Chamaeleon 
I (~ 0.03 pc) and p-Ophiuchus (~ 0.02 pc) to Orion- 
Trapezium (~ 0.002 pc) is largely a projection effect. In 
other words, there appears to be an approximately univer- 
sal distribution of separations in the binary regime, giving 
-/3 bi „ 



N biD 



K blD 9~ 



(6) 



with /3 bin = 2.0 ±0.2 and K bia satisfying the normalization 
condition 



N bi „(9)2n9d9 ~ 1 , 



(7) 



i.e. an average of about one binary (or multiple) companion 
per star. Therefore in denser clusters the intercept between 
the binary regime and the clustering regime shifts to smaller 
projected separations. 

However, the distribution of separations in the cluster- 
ing regime appears to be much less universal. Specifically, 
/3 C ] U , although always well defined, ranges from 0.0 (Simon's 
minimum for Orion- Trapezium) to 0.9 (our maximum for 
Taurus) - see Table 2. Our experiments confirm the con- 
clusions of Bate et al. (1998). The form of N(9) is very 
sensitive to the treatment of boundary effects, and to the 
sample used; this is well illustrated by the work of Naka- 
jima et al. (1998). One possible reason for the very small 
/3 c iu value we obtain for Chamaeleon I (as compared with 
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Taurus) is that we have used an X-ray selected sample here 
and this may favour the diaspora of weak-line T Tauri stars 
(WTTS), as against the more centrally clustered classical T 
Tauri stars (CTTS). By contrast, in Taurus most of the stars 
we use were initially detected in emission-line surveys, and 
therefore we might expect CTTS to be better represented 
than WTTS, giving a large /3 c i u , as observed. It is also clear 
that the procedure of fitting the two power laws simultane- 
ously can make a significant difference to the result. When 
we analyzed Simon's Taurus data by fitting two power laws 
separately, we obtained the same fitting parameters as he 
did. 



2.6 Box dimension 

We have also calculated the box-dimensions of these two 
SFRs. To do this we place a uniform Cartesian grid over 
the star field and count the number of grid-cells M O0 c which 
contain stars. This is repeated for several different values of 
the grid-spacing AO. 

Strictly, the box-dimension is given by 



Dbox - l™o\dln[l/A0]j 



(8) 



However, since the clustering with which we are concerned 
has a finite dynamic range, we have adopted the following 
procedure. 

First we determine the range of AO values to be con- 
sidered. To do this we construct a square box of side A# max 
which just contains the cluster under investigation. If the 
cluster contains M to t stars, we replace them with M to t points 
distributed randomly within the boundaries of this box. Next 
we find the smallest grid-spacing A# min for which the box- 
dimension of the randomly distributed points equals 2. The 
range of AO values is then taken to be A# min < AO < 
A(9 max . 

Finally we return to the star cluster and evaluate the 
slope of the plot of log[7V occ ] against log[A#] for A6 values 
in this range: 



I dln\M,^\ \ 
\dln[\/A0] I 



(9) 



A9 min <Ae<A9 max 



For Chamaeleon I this gives Z) b0 x = 1.51 ±0.12 in the range 
0.29° < A6 < 2.0°; and for Taurus, D box = 1.39 ± 0.01 in 
the range 3° < A6 < 20°. 

Given the small dynamic range of A6, the significance of 
Dbox is simply as a descriptor for the structure on the largest 
scales. It has no significance unless quoted along with the 
dynamic range of A9. In particular it does not constitute 
admissible evidence of fractality. 



3 SURFACE-DENSITY MAPPING 
ALGORITHMS 

Gomez et al. (1993) have generated surface-density maps 
of the young stars in Taurus by applying a kernel smooth- 
ing technique. With this technique they demonstrate the 
existence of several sub-groups of Pre-Main-Sequence stars 
within the Taurus cluster. Such maps are useful for compar- 
ing the distribution of young stars with the distribution of 



dense gas from which stars form, and also for comparing with 
maps generated by numerical simulations of star-formation. 
However, any such map depends critically on the choice of 
smoothing length. It is unclear how Gomez et al. made this 
choice. We present here two algorithms for determining the 
smoothing length objectively. Algorithm I (SCATTER) uses 
a universal smoothing length. Algorithm II (GATHER) uses 
a local smoothing length. We stress that although the final 
maps are evidently similar morphologically, this does not 
mean that the algorithm used to generate a map is unim- 
portant. The information extracted from the raw star field 
is substantially different. The principal differences lie in the 
contour levels, which have a much larger dynamic range for 
Algorithm II (GATHER). 

In both cases a grid is drawn over the star field, and 
the surface-density of stars is calculated at each grid point. 
Standard software is then used to plot contours of constant 
surface-density. Provided the grid is sufficiently fine, the pre- 
cise grid spacing has no significant influence on the final 
map. 



4 ALGORITHM I: SCATTER 

Given A/tot stars with positions (x n , y n ) (n = 1 to A/tot), Al- 
gorithm I evaluates the arithmetic mean (s) of their angular 
separations: 



E E 

n — 1 n' — (n+1) 



Mot (Me - 1) 

| [(x n - x n ,) 2 + (y n - y n ,f] 5 | . 
The universal smoothing length h is then given by 



h = 



1/2 



(10) 



(11) 



.Mot . 

Each star n is considered in turn and distributed to 
the surrounding grid points (x,y) using a normalized linear 
smoothing kernel with compact support, 

3(h - r)/nh a , r < h 

r > h 



W(r) 



(12) 



where r 2 = (x - x n ) 2 + (y - y n ) 2 . 

The contributions from individual stars are summed 
at each grid point to obtain a surface-density array, and 
then contour diagrams are plotted. The minimum contour 
is drawn at 

A mi „ = ^; (13) 

so that isolated single stars are not contoured, i.e. Af min > 
W(0). 

Since the kernel is normalized, the final surface-density 
map is also normalized, in the sense that 



N(x,y)dxdy = M tc 



(14) 



However, by using a universal smoothing-length h, we 
are inevitably capturing only structures on scales of order h. 
Structures on scales much smaller than h are smoothed out, 
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Chaniaeleon I 



o 



-2 




2 -1 

log (e) 

Taurus 



CD 

o 




-2 



Figure 1. Logarithmic plots of N(8) versus 8 for Chamaeleon I (left) and Taurus (right). The crosses with lcr error bars are the corrected 
points which have been fitted with Eqn. (^) (thin solid curve). The open circles represent separations which are ignored in making the 
power-law fits (see discussion in Section 2.3). The filled circles are the uncorrected points. The thick solid lines marked DM show N(0) 
for companions to field stars, as obtained by Duquennoy & Mayor (1991). 
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Table 2. Parameters of the power-law fits to N(8) 



SFR Number binary regime break clustering regime Referena 

of stars 







° L per cleg" J 


/3 bi „ 


Sclb/pC 


i r iCi„ i 


/3clu 




Trapezium 


35 


-2.30 ± 2.52 


2.13 ± 0.65 


0.002 


4.65 ± 0.51 


0.20 ± 0.21 


S97 


Trapezium 


744 




2.00 


0.003 




0.11 ± 0.02 


BCM98§ 


Lupus 


65 


-2.47 ± 0.46 


2.2 ± 0.4 


0.0053 


1.23 ± 0.23 


0.82 ± 0.13 


NTHN98 


Cham. I 


103 


-2.11 ± 0.19 


2.1 ± 0.2 


0.013 


1.36 ± 0.03 


0.57 ± 0.04 


NTHN98 


p Oph. 


96 


-3.42 ± 0.32 


2.5 ±0.3 


0.026 


0.91 ± 0.03 


0.36 ± 0.06 


NTHN98 


p Oph. 


51 


-1.20 ± 0.41 


1.88 ±0.10 


0.024 


1.44 ± 0.24 


0.50 ± 0.19 


S97 


Taurus 





-2.19 


2.15 


0.04 


0.53 


0.62 


L95 


Taurus 


80 


-1.65 ± 0.25 


2.01 ± 0.06 


0.06 


0.54 ± 0.15 


0.64 ± 0.19 


S97 


Cham. 


94 


-2.80 ± 0.16 


2.4 ± 0.5 


0.083 


-0.09 ± 0.02 


0.55 ± 0.03 


NTHN98 


Orion OB 


361 


-1.66 ± 0.13 


1.6 ± 0.4 


0.13 


0.95 ± 0.005 


0.15 ± 0.02 


NTHN98 


Cham. I 
Taurus 


137 
216 


-1.33 ± 0.21 
-1.62 ± 0.10 


1.97 ± 0.07 
2.02 ± 0.04 


0.027 
0.032 


2.00 ± 0.04 
0.54 ± 0.01 


0.28 ± 0.06 
0.87 ± 0.01 


This work 
This work 



a Refs: L95: Larson (1995), S97: Simon (1997), BCM98: Bate, Clarke & McCaughran (1998), NTHN98: Nakajimaet al. (1998). 
b Only results using the McCaughrean et al. (1996) survey are quoted. If boundaries are taken into account, /3 c i u becomes 
-0.02 ± 0.01. 

c Larson used results from several different surveys, but analyzed them separately, so there can be no unambiguous entry for 
his work in the second column of the tabic. 



because the constituent stars are closer together than their 
individual smoothing- lengths. Structures on scales much 
greater than h are not contoured because the kernels of their 
constituent stars do not overlap. If clustering is hierarchi- 
cal (e.g. Larson 1995), the limited dynamical range of the 
SCATTER algorithm will suppress this fact. 

Nonetheless, the resulting contour maps appear to re- 
produce the features which the human eye sees. This is 
demonstrated in Figures 2 and 3. The raw star fields are 
shown in Figure 2, with Chamaeleon I at the top, and Taurus 
at the bottom. The contour maps obtained with the SCAT- 
TER algorithm are shown in Figure 3. Note that the linear 
scale of the Chamaeleon I map (~ 1°) is about ten times 
smaller than the linear scale of the Taurus map (~ 10°). 



5 ALGORITHM II: GATHER 

In Algorithm II we compute for each grid point, the distance 
dk to the fc th nearest-neighbour star, for k < fc max . The 
surface-density N s at that grid point is then calculated by 
gathering contributions from these stars according to 



E 

V fc=l 



W k IT d\ 



k - 



Thus the effective smoothing length is an average of the 
distances d k to the nearest fc max stars. The best results are 
obtained with a weighted sum of the squared distances. 

By invoking an adaptive local smoothing length, this 
algorithm ensures that the computed surface-density is finite 
everywhere - albeit very low in places where the stars are 
few and far between. This extends the dynamic range of the 
algorithm, and thereby increases its chance of capturing any 
hierarchical structure in the underlying star field. 



Algorithm II is objective provided that the choice of 
fcmax and the scheme for computing the weighting factors 
u>k are objective. For simplicity we use: 



w k oc exp [— |fc - 4|] 
with 

fe = fcmax 



(15) 



(16) 



and fc max = 7. Two considerations inform this choice. First, 
very few grid points will have more than 7 natural-neighbour 
starsp]. Second, we do not wish to weight heavily the closest 
neighbouring stars, otherwise we will over-emphasize binary 
systems. 

However, the resulting surface-density map is only nor- 
malized, in the sense of Eqn. (Q, if the surface-density is 
truncated at a finite value. For Chamaeleon I this truncation 
value is 7V min = 47.7 stars per square degree, and for Taurus 
it is iV min = 4.06 stars per square degree. In other words, at 
all grid points where N g falls below the truncation value, we 
set N g to zero. 

Again the resulting contour maps appear to reproduce 
the features which the human eye sees. This is demonstrated 
in Figure 4. Inspection of the limiting contour values listed 
in Table 3 shows that the GATHER algorithm does indeed 
have a significantly larger dynamic range than the SCAT- 
TER algorithm. For this reason GATHER is our prefered 
algorithm. 



^ If the grid point and its nearby stars are treated as nodes, then 
the natural-neighbour stars of the grid point are those stars whose 
Voronoi cells share a boundary with the Voronoi cell of the grid 
point, e.g. Braun & Sambridge (1995) 
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Figure 2. The star positions for Chamaeleon I (top) and Taurus (bottom). The plots are centred at a = 11 hrs 7 mins 45.6 sees and 5 
= - 77° 5' 24" for Chamaeleon I and at a = 4 hrs 27 mins 2.4 sees and <5 = 24° 33' for Taurus. 
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Figure 3. The surface-density plots generated by SCATTER, for Chamaeleon I (top) and Taurus (bottom). The plots are centred at 
= 11 hrs 7 mins 45.6 sees and <5 = - 77° 5' 24" for Chamaeleon I and at a = 4 hrs 27 mins 2.4 sees and <5 = 24° 33' for Taurus. 
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Figure 4. The surface-density plots generated by GATHER, for Chamaeleon I (top) and Taurus (bottom). The plots are centred at a 
= 11 hrs 7 mins 45.6 sees and <5 = - 77° 5' 24" for Chamaeleon I and at a = 4 hrs 27 mins 2.4 sees and <5 = 24° 33' for Taurus. 
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6 RESULTS 

Contour maps obtained by the two algorithms described 
above are shown in Figures 3 and 4. On all four maps, 
the contours are equally spaced in logarithm between iV min 
and A?p 0ak , but since A pcak corresponds to the highest sin- 
gle pixel, the maximum contour drawn falls below -/V poak at 
-/V max . On all maps, only 5 contours are drawn, to avoid con- 
fusion. 

We see that Chamaeleon I divides into two main sub- 
clusters, having > 100 stars per square degree, and linear 
extent ~ 0.4° or ~ 1 pc. Within these are three more com- 
pact sub-sub-clusters having > 500 stars per square degree, 
and linear extent ~ 0.1° or ~ 0.25 pc. Taurus contains 
many more sub-clusters. Gomez et al. (1993) identify six 
sub-clusters, but of course the precise number depends on 
how a sub-cluster is defined. If we require a sub-cluster to 
be an isolated peak of at least one order of magnitude in 
surface density, we find 7 sub-clusters in the GATHER map 
- which is the more discriminating map. The sub-clusters in 
Taurus are much less dense (typically ~ 10 stars per square 
degree) than those in Chamaeleon I. 

In particular we note that the two Chamaeleon I sub- 
clusters are similar in linear size and separation to the two 
sub-clusters in the north-west corner of Taurus (Groups I 
and II of Gomez et al. 1993, at a ~ 4hrs 15mins and 8 ~ 
28°). However, the two sub-clusters in Chamaeleon I contain 
about five times as many stars as Groups I and II in Taurus. 

In fact the sub-clusters in Chamaeleon I are probably 
bound. Each sub-cluster contains about 50 stars, so its total 
mass is ~ 10 g, and each sub-cluster has a radius ~ 0.5 pc. 
Therefore they are bound provided their intrinsic velocity 
dispersion is Av 1 km s _1 . 



7 CONCLUSIONS 

We have calculated the surface-density of companions per 
star as a function of angular separation for Chamaeleon 
I and Taurus. We have used a fitting procedure which 
avoids the subjectivity of locating the elbow by eye. For 
Chamaeleon I , we obtain N(9) oc 9~ 197 for 9 <C 0.013° 
and N(9) cx 0" - 28 for 9 > 0.013°. For Taurus, we ob- 
tain N{9) oc <r 2 ' 02 for 9 « 0.011° and N{9) oc (T - 87 for 
9 > 0.011°. 

In combination with the results of Larson (1995), Si- 
mon (1997) and Nakajima et al. (1998), our results suggest 
that universal processes are at work determining the mul- 
tiplicity of newly formed stars, with f3 bin ~ 2.0 ± 0.1, and 
K hin ~ 3.±l.x 10" 2 stars per square degree. One possibility 
is that the statistics of binary and multiple systems are de- 
termined by dynamical fragmentation occuring in gas which 
has suddenly become dense enough to be cooled efficiently 
by dust (Whitworth, Boffin & Francis 1998). 

In contrast, the large-scale clustering of young stars 
shows quite a large variation from one star-formation re- 
gion to another. In general the power-law fit to N(9) in the 
clustering regime is rather tight, and hence the values of 
and /3 du obtained are well determined, although they 
do depend quite sensitively both on the sample used, and on 
the analysis technique used. The differences between star- 
formation regions are not solely due to different overall sur- 



face densities. There is no compelling evidence for a univer- 
sal clustering process at work. 

We have also calculated the box-dimensions for the two 
star clusters: 1.51 ±0.12 in the range 0.29° < A9 < 2.0° for 
Chamaeleon I, and 1.39 ± 0.01 in the range 3° < A9 < 20° 
for Taurus. We emphasize that these are only descriptors of 
the clustering and not admissible evidence for fractality. 

We have presented two objective algorithms for con- 
structing surface-density maps of star clusters. The SCAT- 
TER algorithm has the advantage that it is implicitly nor- 
malized, but it has limited dynamic range and therefore sup- 
presses hierarchical clustering. The GATHER algorithm has 
the advantage of increased dynamic range, but it has to 
be normalized explicitly. Both algorithms capture the struc- 
tures seen by the human eye, and both avoid contouring 
individual stars 

The surface-density maps show that the sub-clusters in 
Chamaeleon I are very similar in linear size (~ 1 pc) and sep- 
aration (~ 2 pc) to those in Taurus, but apparently much 
denser; if we assume that the sub-clusters in both SFRs 
are approximately spherical, then those in Chamaeleon are 
denser than those in Taurus by a factor i> 5. There ap- 
pears to be a characteristic length-scale for cluster formation 
which is roughly independent of cluster mass. For instance, 
Lada et al. (1991) find that the embedded clusters in L1630 
also have sizes and separations of this order. 

One of the interesting problems which Bate et al. (1998) 
identify is that a given /3 du could derive from a fractal distri- 
bution, but it could equally derive from single-level cluster- 
ing. The resolution of this dichotomy may lie in the analysis 
of surface-density maps of the type we have generated here. 
For instance, in principle fractality should be revealed by a 
plot of perimeter versus area for the iso-surface-density con- 
tours. However, in practice, we have found that in Taurus 
and Chamaeleon I there are insufficient contours to obtain 
a reliable dimension. And in reality, the notion of fractal 
dimension may not be very meaningful for a property (the 
surface density) with such a small dynamic range, or a star 
field with so few stars (~ 100). 

Another problem which Bate et al. (1998) identify is 
that N(9) does not give any indication of the density pro- 
files of individual sub-clusters. Here surface-density contours 
can clearly help, since even small- A clusters are resolved in 
unconfused regions like Taurus. For an isolated cluster de- 
fined by contours at surface-densities N k (with N k+ i > N k ) 
and enclosing areas A k , we have 



dln[N] 



dln[R] 



lim 

fc — >fc max 



Aj /2 (N k - Afc-i) 



and 



dln\n] 



dln[R] 



dln[N] 



dln[R] 



+ 1 



(17) 



(18) 



where fc max is the central contour in the cluster, N is the 
surface-density, and n is the volume-density. 
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Table 3. Contour limits and peak values for the maps in Fig. 2 

Algorithm Chamaeleon I Taurus 





JV mm 


range 


-/Vmax 


AW 


JV min 


range 


-/Vmax 


JVpeak 


SCATTER 


140 


x8.2 


751 


1142 


3 


xll.7 


21 


35 


GATHER 


47.7 


x287 


4413 


13690 


4.06 


x256 


342 


1038 



On all plots there are 5 contours, and the contours are equally spaced logarithmically between N, 

and JVpeak- 



94157. We thank the referee for suggesting that we make Whitworth A. P., Boffin H. M. J., Francis N., 1998, MNRAS, 299, 

an explicit comparison between our distributions and the 554 

results of Duquennoy & Mayor for field stars. Whitworth A. P., Chapman S. J., Bhattal A. S., Disney M. J., 

Pongracic H., Turner J. A., 1995, MNRAS, 277, 727 
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